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The renormalization of the magnetic field penetration depth A owing to the electron-electron 
correlation is discussed with its application to high-T c cuprates. The formula for the current 
carried by quasiparticle with the Umklapp scattering is derived, on the basis of which we 
investigate how the value of A -2 deviates from that of n/m* where n and m* are the carrier 
density and the effective mass respectively. Although this deviation is small in the case of weak 
momentum dependence of the vertex, this is large and negative owing to the non-negligible 
value of the backflow in the case of the strong antiferromagnetic spin fluctuation. The observed 
doping dependence of A~ 2 in high-T c cuprates, specifically a peak structure at the slightly 
overdoped region, is explained by the analytical consideration and the numerical calculation 
based on the perturbation theory and the spin fluctuation theory. The consistency between 
A~ 2 and dA _2 /dT at absolute zero, which is the problem the isotropic model fails to explain, 
is also obtained by our theory. 
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§1. Introduction 

In high-T c cuprates some transport phenomena show characteristic behaviors, especially in the 
optimal and the underdoped regions which do not seem to be explained by the conventional Fermi 
liquid theory. For examples, the resistivity is proportional to the temperature T, rather than T 2 .0) 
The Hall coefficient is temperature dependent and increases as the doping level decreases.^ It 
is also seen that not only these transport phenomena, but also the magnetic properties like the 
nuclear spin-lattice relaxation rate and the one-particle spectrum show the pseudo gap phenomenon 
which is the gap like behavior with no long range order. The magnetic field penetration depth A 
in the high-T c cuprates attracts attention because A is long in the underdoped region where the 
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pseudogap phenomena appear and then some authors suggest the low carrier density based on the 
relation A~ 2 oc n/m* where n,m* are the carrier density and the effective mass, respectively,!^ and 
a connection between long A and the pseudogap phenomena.!^ 

In the above phenomena it is found in the early stage of the study on the high-T c cuprates that the 
T-linear resistivity is explained by the fact that the lifetime of the quasiparticles is proportional to T 
in the presence of the strong antiferromagnetic spin fluctuation.!)'!^ This means that it is sufficient 
to consider an one-particle property like the imaginary part of the self-energy. The pseudo gap 
phenomena can also be explained by the superconducting fluctuation with the two dimensional 
and the strong coupling nature of the system.B'i' In this case the absolute value of the imaginary 
part of the self-energy shows the maximum at the Fermi level and therefore one particle spectrum 
decreases. On the other hand it is found that the interaction between the quasiparticles is essential 
for the explanation of the Hall coefficient M This work clarifies the importance of the interaction 
between the quasiparticles besides one-particle properties like the self-energy. 

As for the magnetic field penetration depth in high-T c cuprates, there has been no satisfactory 
theories until now. As noted above it seems that A~ 2 is proportional to the doping concentration 5, 
and from the combination with A~ 2 oc n/m* it seems to be reasonable to expect that the effective 
carrier density is identified with the doping concentration (i.e. n oc 5)0 The theories based on the 
t-J model naturally accept this conjecture because of no double occupation. There is no need to 
explain why A is so small but it needs only to take it as an external parameter in these theories. 
However the doping dependences of A -2 show that A -2 increases with 5 in the underdoped and 
optimal regions, and it has a peak at the slightly overdoped region, and then it decreases as 5 
increases in the overdoped regionB It is obvious that the latter behavior cannot be explained by 
the theory based on the t-J model. On the other hand n should be proportional to 1 — 5 in the Fermi 
gas theory because 5 = means the half-filled, and therefore based on A~ 2 oc n/m* it can explain 
the overdoped region but cannot explain the optimal and underdoped regions. If we take account 
of the interaction between the quasiparticles, it is known that A~ 2 oc nil + Ff/3)/m*^ Here Ff 
is one of the Landau parameters and this expression is based on the isotropic model. Though the 
behavior of this parameter is not known, it is seemingly possible to explain the doping dependence 
of A~ 2 by controlling the parameter Ff. However it is shown that this parameter theory fails if we 
try to maintain the consistency between the values of A~ 2 and dA _2 /dT at absolute zeroJlil These 
considerations indicate that the isotropic formalism presented by Leggett should be reexamined to 
investigate the realistic metallic system like the high-T c cuprates, and it cannot be permitted to take 
Ff as an external parameter either, necessitating the investigation of the quantities corresponding 
to Ff in detail. 

The Umklapp scattering, a main interest in this paper, is important for the transport phenom- 
ena. In the previous paper a systematic discussion on how the many body effect appears in the 
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superconducting state, is presented S3 Some examples presented there, are the magnetic field pen- 
etration depth and the optical conductivity. At absolute zero the expression of the magnetic field 
penetration depth contains the interaction between quasiparticles in the same way as in the cy- 
clotron resonance frequency© and the Drude weight ©'0 In ref.18 it is shown that the absence 
of the Umklapp scattering leads to the infinite conductivity. Although this fact can be shown by 
using the Ward-Takahashi identity related with the momentum conservation, the essential point 
in ref.18 is that the formula for the resistivity is written explicitly by the Umklapp process. The 
formula for the current which is explicitly written by the Umklapp scattering has not been known. 
The derivation of this formula is given in this paper. With regard to high-T c cuprates it is needed 
to clarify the relation between the Umklapp scattering and the antiferromagnetic spin fluctuation 
and it is shown that the strong antiferromagnetic spin fluctuation leads to the strong Umklapp 
scattering on the current. 

In §2 the formula for the current carried by quasiparticles written by the Umklapp scattering 
is derived and the other property of the current is briefly discussed. In §3 the model and the 
approximations for the analytical and the numerical calculations are presented. Some discussions 
on the influence of the superconducting transition is also given and the validity of the calculation 
of the current in the normal state is shown. In §4 the analytical investigation for the current 
is presented in the two specific cases. One is the case where the momentum dependence of the 
irreducible four point vertex is weak and the other is the case with the strong antiferromagnetic 
spin fluctuation. It is shown that the discussions based on A~ 2 oc n/m* is approximately valid 
only in the former case. In the latter case the behaviors of the current differ with the others 
depending on the points on the Fermi surface and it is shown that the current on the points with 
the strong antiferromagnetic spin fluctuation is decreased much. How the effect of the Umklapp 
scattering on the current depends on the irreducible four point vertex is also presented and it is 
shown that it is largest in the case of the strong antiferromagnetic spin fluctuation and is smallest 
in the case of the ferromagnetic spin fluctuation. It is also presented that the backflow in the case 
of the strong antiferromagnetic spin fluctuation has an opposite sign to the bare velocity. In §5 the 
numerical calculations with the self-consistent second order perturbation theory (SC-SOPT) and 
the fluctuation exchange approximation (FLEX) are presented. The results which are consistent 
with the experiments in the respects of the doping dependence of A -2 at absolute zero and the 
consistency between A -2 and dA~ 2 /dT, are obtained. In §6 a brief consideration on the relation 
between the superconducting fluctuation and A -2 is given and it is indicated that the small A~ 2 
at absolute zero does not directly mean the large thermal fluctuation. In §7 the summary and the 
discussion are given. 

In this paper we put fi = k-Q = c = 1. 
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§2. Current Carried by Quasiparticles 

The current carried by quasiparticles plays an important role in physical quantities of the colli- 
sionless region in the normal state, for example, the cyclotron resonance frequency and the Drude 
weight. This is usually written by using only quantities defined on the Fermi surface as 

^ = 4, + ^E4,fe'^H v (2.i) 

k 

Here /i is the index of spatial dimensions, V is the volume of the system, £jf is the dispersion of the 
quasiparticles including the chemical potential, vt is the renormalized velocity of quasiparticles 

( v % = ~dle^ an< ^ fk k' = z k^% k' z k' * s ^ e interaction between quasiparticles (F^ is the cj-limit of 
reducible four-point vertex and the notation about w-limit and /c-limit is the same as in ref.19). In 
the superconducting state the magnetic field penetration depth is also written by using this quantity 
at absolute zero. As derived in the previous paper, the magnetic field penetration depth A at 
finite temperature is written as 

Here J FS dSk is the integral over the Fermi surface and Y(k;T) is Yosida function 

nk-,T) = jdtl(- 9 ^y (2.3) 

f(x) := l/(e X//T + 1), = + AjL), ^k ls tne superconducting gap and vl(k; T) satisfies the 

following integral equation, 

*5(fc; t) = m - / Fg J^ {n f k ,k' Y ^ T ^ ^ 

or 

v* u (k;T) = vt(k) + 47r ^ (fc0| ^fc,fc'(l ~ Y(k';T))v:(k';T). (2.5) 

Then it can be seen that v*(k;T) = at absolute zero because of Y(k;T = 0) = 0. Here 

A k k' = z k^fo k' Z k' ^k k' lS ^ e ^" nm ^ OI " ^ ne reducible four point vertex) and this satisfies the 
following equation: 

4,fc' = A kM + yL A k,y y — -^r- J fk",k'- ( 2 - 6 ) 

Here we briefly explain why A" 2 is not generally written as n/m*, and how the many body effect 
enters. The electromagnetic response kernel is written as 

K, v = -J k V k ^GG + FF) k [e)K kv {e)-(^^ (2.7) 

Here G and F are the normal and the anomalous Green's function, respectively, in the superconduct- 
ing state in the usual sense, 19 ), m is the bare mass, v k is the bare velocity and Ay^(e) is the three 
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point vertex which satisfies the integral equation with the interaction included. J fc := J J2k ^ s 
used in the zero temperature formalism hereafter. In the case of the finite temperature formalism 
Ik := / 1^7 Sfc is used. The first term of the r.h.s. of eq.(2.7) is often called paramagnetic term 
and the second term is called diamagnetic term. In the usual textbook© it is explained that the 
paramagnetic term vanishes at absolute zero in the superconducting state. This is valid in the case 
where the superconducting gap grows in the whole Fermi sea and then the integral which includes 
GG+FF reduced to Yosida function and vanishes at absolute zero. However in all superconductors 
the superconducting gap grows only near the Fermi surface and the incoherent part contributes 
to the integral which does not vanish in the lattice system. Then the eq.(2.7) is reduced to the 
following equation at absolute zero, instead of to —n/m. 

K,u = -j k V k ,{GGY£{e)k kv (e) - . (2.8) 

In the same way A^. (e) satisfies the following integral equation. 

Afa,W =v k „ + ^J fcfV (6,0(GG)Jf(0A^(0. (2-9) 

(1^ £.'( e ' e ') is the irreducible four-point vertex.) Then from eqs.(2.7) and (2.9) the formalism of 
magnetic field penetration depth based on Fermi liquid theory is derived by integrating the high 
energy part to derive the low energy dynamics as investigated in detail in ref.14. In the specific case 
where the momentum is conserved it can be shown that the paramagnetic term exactly vanishes 
at absolute zero as follows. In this case Ward-Takahashi identity which reflect the momentum 
conservation is written as 

/ 0£?(e)\ 

Afa,(*) = " -J^J (2-10) 

(S'^(e) is the normal self-energy.) This relation holds at arbitrary value of e. Then the paramagnetic 
term is 

= 0. (2.12) 

This is the direct proof of the disappearance of the paramagnetic term in the case where the 
momentum conservation holds, which is not via the Fermi liquid parameter like Ff. 

Next we discuss the role of the Umklapp scattering. Hereafter in this section the discussion is 
given by using the Green's function in the normal state. In the superconducting state equations 
become more complicated, but the main results are not altered as partly discussed in ref.14 about 
the renormalization factor and the renormalized velocity and in §3.1. By using the fact that the 
renormalized velocity and the renormalization factor is written as (k is a simplified notation of 
(fc, e) in the zero temperature formalism), 
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and 

z k 1 = l + J k Tl k ,{^e')[G{k'fr (2.14) 
respectively, and also noting that 

[G(kfr -\G{kf} k = 2mS(e)S(e k ), (2.15) 

eq.(2.1) is written as 

fk, = «*„ + *k I r fc,fc'(°' OP( - « fc/i ), (2.16) 

or 

Here d ^jp is the same notation as in ref.21 and it means the derivative with the Fermi surface also 
transformed, in the contrary to 9 ^ k ^ , the derivative with the Fermi surface fixed. Reflecting the 
simultaneous transformation of the Fermi surface, corresponds to the w-limit, and to 

the fc-limit owing to the inclusion of the discontinuity on the Fermi surface. 
The quantities, d ^j^ and 9 , satisfy the following integral equations, respectively. 

dS(fc ) / ^^^^ ^. , d£(fc') 
dfc 

and 



//(^^)[^') 2 r(^ + ^P) (2-18) 



9m ~ ~ I I(k, k')[G{k'fr ( 1 - ^r 1 ) ■ (2-19) 



de Jk' V de' 

dS(fc) dE(fc) , 



Then it is verified that w k (e) := ; H — ^~ v k satisfies the following integral equation. 

= «fc(<0 + / ^[G^OT™^') (2.20) 

Jk' 

= u k (e)+ f T»{k,k>)[G{k>?ru k ,{e'). (2.21) 

Jk 1 

Here 

«*„(«) = I k')[G(k'fr (l - ^f-) (v^ - v klt ). (2.22) 

This is rewritten as 

= X i |r(t ' k ''' ,)|2<3(t " ,)G(F + ?) (~^f^) + " " <2 ' 24) 

This equation shows that «^.(e) = o in the absence of the Umklapp process and the concrete 
examples are given in §3 in the case of SC-SOPT and FLEX approximations. By using the above 
quantities the backflow term is written by the difference between the two kinds of the momentum 
derivative of self-energy. 

am_dm_B k 



dk dk 



z k 
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Here is the backflow term. 



In ref.15 is calculated based on second order perturbation theory with respect to the on-site 
Coulomb interaction. It is mentioned there that the four-point interaction vertex is taken to satisfy 
the conservation law and they obtained the results that the vertex correction make the positive 
value in some parameter regions. If this holds, it follows that the value of current carried by 
quasiparticles exceeds the value of the bare electron's velocity. However it is not considered to be 
possible for the interacting electrons to carry the current which exceeds the non-interacting case 
because if it were possible, gathering the momentum from the crystalline lattice and carrying it 
would occur. If we see the general expressions above it is obvious that their treatment is incorrect. 
To satisfy the conservation law and get the correct value of , it is needed to use eq.(2.24) as the 
vertex correction with \T(k, k'\ q)\ 2 replaced by U 2 in the case of the second order perturbation. 
The fact that j-^ does not exceeds t>^ is made clear by eq.(2.16) with the notion that is an 
odd function with respect to k and f k , ^'(0, e')[G{k') 2 \ JJ > in the Fermi liquid system. Then 
j'k < Vfc^ in k with Vf.^ > 0, and this is the basic inequality which can be obtained from the 
above general discussion. This inequality also guarantee that A -2 is always smaller than 47re 2 n/m. 

By taking the above results into account we consider how the currents j* in the lattice system 
changes from the value j* = v which is universal in the case of the continuum. If we consider 
this correction in terms of the renormalization factor (z), it can be seen that this quantity is 
changes as (1 — z~~ 1 )/z~ 1 plus the extra contribution proportional to z which depends on the 
system considered. The presence of this extra contribution violates the relation A oc n/m* and 
expresses one of the characteristics of the system like the momentum dependence of the interaction 
between the quasiparticles. This property in high-T c cuprates is investigated in detail in §4 and 5. 

§3. Model and Approximation 

We take the following Hubbard Hamiltonian with on-site Coulomb interaction. 



k,a k,k',q 
Here £^ is the dispersion of electrons which is taken similar to that in ref.9. 

£^ = —2t(cosk x + cosky) + 4t' cosk x cosk y — 2t"(cos2k x + cos2k y ). 




(3.1) 



(3.2) 



In the numerical calculation we put t'/t 
The self-energy within SC-SOPT is 



0.16, t"/t = 0.20, U/t = 7.0 and t = 1. 




(3.3) 



Here k = (k, ie n ) (e n is the Matsubara frequency), 



G{k) = 7 



1 



(3.4) 
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and 

x(q) = -^J2 G ( k ') G ( k ' + iy ( 3 - 5 ) 



V V 



The irreducible four point vertex is 



I(k,k') 



SG(kf) 

= 2 X (k-k') + (/)(k + k'). (3.6) 
The first line of this equation follows the conserving approximation by Baym and Kadanoff Here 

^) = -^E G ( fc ')G(g-A:'). (3-7) 

k' 

r^,,(e, e') and T^,,(e,e') in the basic formalism in Appendix are given by the following expres- 

lt> {e,e') = 2 X l h ,(e-e>) (3.8) 



k,k' 
sions 



and 

T k^ e ' )= ^ + k^ + ^ ( 3 - 9 ) 
By using the above I(k,k'), Wfe(e) which expresses the Umklapp term is explicitly written as 

u k& = U 2 j k , j q G k-q(t ~ ")G k > +q {t + u)—l^(v k + v k , - v k , +q - v k q ) (3.10) 

As can be easily seen by explicit calculation —v kq term and ~ v k ' + g i n the above equation 
originate from 2x(q) term in I(k,k'), term from (p(q) and v k term from dTi{k)/de. 
The basic equations in the FLEX approximation are (for an example, see ref.23), 

s (*0 = £E^(?) Gf ( fe -9)> ( 3 - n ) 
The irreducible four-point vertex in this approximation is 



I(k,k') 



6G(k' 

V ( k ~ k ') ~ | E ^[G^ - - g) + G(tf + g)G(fc - ?)]. (3.13) 
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Here 

^H^^^+HttW" 172 (3 - 14) 

r^,,(e, e') and T^^, (e,e') in the basic formalism in Appendix are given by the following expres- 
+ tanh^G^ q (e - W )ImG*, _ q {e> - u)] (3.15) 



sions 



S 



and 

+ tanh^Gj^e - u )ImG*, +9 ( £ ' + w)]. (3.16) 
By using the above I(k, k'), Wj^(e) which expresses the Umklapp term is explicitly written as 

Wfc(£) = L I W ^ )G k-q( € ~ ^ G k' + q^ + ^)^7^K + «fc' " v k>+q ~ V k q) (3-17) 

As can be easily seen by explicit calculation —Vfc_q term in the above equation originate from 
V(q) term in I(k, k'), ~ v y + q term from the second term of r.h.s. of eq.(3.13), Vj^> term from the 
third term of r.h.s. of eq.(3.13) and term from <9£(A;)/<9e. 

By using the above approximations, it is shown that the momentum dependence of Vq(0) is 
remarkable at q = (vr,7r). On the other hand, although Xg(0) has a small peak around q = (vr, 7r) 
and 4>^(0) atq= (0,0), it can be seen that the momentum dependence of the irreducible four point 
vertex in SC-SOPT is rather weak as compared with that in FLEX. 

3. 1 Influence of Superconducting Transition 

Here we briefly discuss how the electronic properties differ between the normal state and the 
superconducting state, and take the renormalization factor and the damping rate of the quasipar- 
ticles as the representative quantities because the renormalization factor is related to the vertex 
corrections to the current by the real part of the irreducible four point vertex and the damping 
rate of the quasiparticles is related to the vertex correction in the hydrodynamic region by the 
imaginary part. It is known that the imaginary part of the vertex caused by the electron-electron 
interaction becomes small at the low energy in the superconducting state because of the decrease 
of the scattering accompanied by the gap formation, while the real part is not so affected at low 
energy except for the uniform magnetic susceptibility in the case of the singlet pairing. 

The imaginary part of the self-energy and the renormalization factor calculated numerically in 
the FLEX approximation is shown in Fig. 1. From this figure it is shown that — ImS^(O) rapidly 
decreases below T c as expected from the above consideration. A notable point concerning the 
damping rate, which is not related to the present discussion, is that the damping rate at the cold 
spot decreases linearly with the temperature and that at the hot spot decreases as y/T. The former 
point causes the T-linear resistivity and the latter point does not affect the conductivity because 
the point with the smallest damping rate mainly contributes to the conductivity. 

On the other hand the renormalization factor almost remains constant as expected also, while 
the change in the spectrum of the spin fluctuation makes the minimum of this quantity at T = T c . 
These considerations justify the calculation of the current in the normal state to estimate the 
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Fig. 1. — ImE^(O) and the renormalization factor (2) at the cold spot (cs) and the hot spot (hs). 5 = 0.10 and 
T c = 0.024. 



magnetic field penetration depth because the vertex correction for the current is the type of the 
function analogous to the renormalization factor. 

§4. Analysis of Umklapp Term 

4-1 The case where the momentum dependence °f Ifcfc' is weak 

This case occurs when the shape of the Fermi surface is not so peculiar and the perturbation 
scheme is valid. Then it is easy to derive the following relation by using the fact that the velocity 
is an odd function of k. 



w k( e )- u k( e ) 



Of 



>k- 



Then the current is written as 



J k, 



- V ka + Z k 



Z k V k, 



de 1 



-o v k, 



(4.1) 
(4.2) 

(4.3) 

(4.4) 
(4.5) 



This equation shows that the Umklapp term reduces to (zu — 1)^^- The current and A -2 is only 
reduced by zu and then the latter quantity is interpreted by writing n/m* as used apriori in the 
some papers Ji^ where m* is the thermal mass which is enhanced by 1/z. However it is considered 
that in the case of high-T c the validity of this perturbation scheme is restricted in the overdoped 
region. This assertion is understood by noting the behavior of the nuclear spin relaxation rate and 



10 



the resistivity measurements. Therefore the use of the form n/m* in the underdoped cuprates, 
which is taken in some papers, is not warranted. 



4-2 The case of strong antiferromagnetic spin fluctuation 

This is the case which is relevant to the optimal and underdoped cuprates. In this case Wq(uj) has 
a sharp peak at q ~ Q and u ~ (Q = (tt, 7t)). The fact that Wq takes a large value at a specific 
momentum reduces the integral including Wq to the one that is the integral of the odd function 
Vfc, and then the term including Wq can be neglected, as confirmed numerically. Then eqs.(2.20) 
and (2.22) (the corresponding equations in Appendix) is reduced to the following equations. 

Hfe(0) = «2(0)+ y^[coth-Im^_ fc ,(- e * +tanh- ±* ImG«(e>)]z k ,w«,(0) 

(4.6) 

and 

«g(0) = y y [cath-Im^C-O ^ + tanh- ^ ImG* (e')](^ - v k ). 

(4.7) 

Here in the former equation we put w?,(e')/(l — dHy{e')/de') as z k 'W?/(0), because only the region 
of small e' contributes to the integral owing to the factor dG(e')/de' . The consistency between the 
numerical calculation and the analytical discussion in this section also verifies this replacement. 
As often taken in the spin fluctuation theories, we use the following approximation to Vq(u>), 

(£ and LO s f are the correlation length and the characteristic frequency of the antiferromagnetic spin 
fluctuation, respectively.) The following form of GjjJ (e) is also used, 

*iRt^\ _ z k 

Then the integration in eq.(4.7) is carried out at T — > and it can be shown that the term which 
includes ImV 1 * vanishes and only the term which includes ImG^ remains. The result is the following 
equation, 

u k [U) - 2 (2vr)2 UF e J FS < i/e + (k-k'~ qy ■ (4 - iUj 

Similarly we obtain 

r> R/ , TTZ kX (Q) f ,,1 W t>(°) 

hew = "gc> + y ra dt < i/p + (fc - g - q? • < 4 - u) 

(vf is the velocity at the Fermi surface.) Then in the limit of £ — > oo the integral in these equations 
reduces to the delta function <5(fc — k' — Q). However, as £ is not infinitely large in real systems, 
we introduce the reducing factor. Then the following equations are obtained. 

w k(°) = + ^ c «™Lq(°)< ( 4 - 12 ) 



Gjfc) * : ,* fe ■ (4-9) 
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and 

«g(0) = c(av k _ Q - v k ). (4.13) 

Here c := 2 (2*)^pg anc ^ a ^ s some constant which is the reducing factor owing to the momentum 
dependence of Vq(0) and < a < 1. By using (v k _Q x , v k _Q y ) = — {v ky , v kx ), the above equations 
are solved as, 

u? (0) - zcauf (0) 

«°>- i -(«,).*' (4I4) 

and 

«2x(0) = -c(«fe e + «i;fcv)- ( 4 - 15 ) 

The equation for y-component is written by exchanging x with y in the above equations. To see 
the behavior of the current on the Fermi surface we analyze these equations on the typical three 
points, at two hot spots (hsl and hs2) and at a cold spot (cs). The locations of these points arc 
shown in Fig.3 of §5. 

«g x (0) ^ -cv kx ,u^ y (0) ~ -cav kx (at hsl), (4.16) 

U L(°) - U ky(°) ~ ~< l + a ) V kx ( at CS )' ( 4 - 17 ) 

u L-- cau fei/' u fc»-- cl, fci/ ( aths2 )- ( 4 - 18 ) 
By using the above equations and z ~ (1 + c) _1 , c >> 1 we obtain 

fk x ^ z k v k x hsl), (4.19) 
i*kx- z k v k x ( atcs )> (4-20) 

From the above equations we can see that the growth of the antiferromagnetic spin fluctuation 
described by a affects mainly near hs2. Then we can see that the reduction of 1/A 2 accompanied 
by the spin fluctuation is mainly caused by hs2. It can be said that this point is strongly affected 
by the Umklapp scattering as can be seen in fig. 2. 

On the other hand the temperature dependence of 1/A 2 at the low temperature is dominated by 
the point cs and dA~ 2 /dT| T=0 is determined by j* 2 unlike j* in the case of A~ 2 at absolute zero. 
These properties are shown in the formula of dA _2 /dT at absolute zero derived from eqs.(2.2) and 
(2.4) as 

In this equation the main contribution comes from near point cs because of the pairing symmetry 
of high-T c cuprates. Then this equation indicates that the reduction of 1/A 2 (0) accompanied by 
the spin fluctuation does not mean that the reduction of (d/dT)A~ 2 (0) because j* at cs is not 
affected by the growth of the spin fluctuation as can be seen from eq.(4.20). This is also confirmed 
numerically in §5.3. 
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4-3 Nature of Umklapp process 

The large value of the Umklapp term in the case of strong antiferromagnetic spin fluctuation is 
understood by considering the extreme example as follows. If Vq(ui) and Wq(uj) are, respectively, 
replaced by VS(q — Q)5(oo) and WS(q — Q)6(u>), then u k (e) is written as 

Uk(e) = V L ^-Q^k'-Q^^r 1 ^ - v k q) 

= W J w G k _ Q (e)G k ,_ Q (e>)—l^v k . (4.23) 

Here v k q is replaced by — v k . 

As the q-dependence of Vq(uj) and Wq(u>) is weakened, it is easily seen that the magnitude of 
the above u^(e) becomes small. On the other hand if Vq(u) and Wq{io) have sharp peaks at q = o 
uj = 0, i.e. in the case of strong ferromagnetic fluctuation, (though this is not the case of high-T c , 
we consider it as just an example to illustrate the effect of the Umklapp scattering) it can be seen 
that the Umklapp term vanishes as below. 

«fcW = y I G k ^G k ,(e>)-l^(v k - v k ) 

= 0. (4.24) 

These behaviors can be understood intuitively by drawing the processes like in Fig. 2. The 
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Fig. 2. The characteristic scattering processes with and without the Umklapp scattering 

scattering process is (k, k') — ► (k — q,k' + q), and then by using the notation in this figure the case 
of the strong antiferromagnetic spin fluctuation corresponds to the case where the process of A 
dominates and the case of the strong ferromagnetic spin fluctuation corresponds to the case where 
the process of B is dominant. The process A is clearly the Umklapp process and the process B 
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means that the electron goes back to the beginning point on the Fermi surface and this process has 
no contribution to Ufc(e). Both processes are equally contained in the case where the momentum 
dependence of the irreducible four point vertex is weak. In the actual system both processes exists 
and which type of these processes is dominant is the factor for determining whether the current 
j* < v* or not. 

4-4 Backflow 

In §2 it is shown that the backflow can be estimated from the two kinds of the derivative of 
self-energy by k, <9£^(0)/<9fc and dS^(0)/dfc. The integral equations which are obeyed by these 
quantity is written as follows; 

d££(0) r e' D dReG*{e') e > OReV* ,,(-e') 

_^i = / [ cot h4-ImK R .,(-e') -J^+tanh^ ^ ImG^e') 

dk Jk> 2T k-k v ' de' 2T de' k v ' 

+ d? ( tanh 2r) ^ v tk'^') l ^ G k'^ [ v k' + -k^J ( 4 - 25 ) 

and 



dfc 



/ [coth— ImF* ,,(-e') +tanh- ImG^e')]^ v y + 

Jk> IT k-k y ' de' 2T de' k y n k \ k 



Here S^,(e') is replaced by S^,(0) by considering that small e' region dominates in the integral as 
before and terms including Wq(u) can be neglected for the same reason in the previous subsection. 
The main difference between the above two equation is the presence of the third term in eq.(4.25) 
and this is the difference between the k and cj-limits in the finite temperature formalism. The 
integral of the first term and the second term in eqs.(4.25) and (4.26) are carried out in the same 
way in eq.(4.7) and it can be seen that the third integration in eq.(4.25) has opposite sign and twice 
in magnitude as compared with the second one. Then the eqs.(4.25) and (4.26) reduce to 

asf(o) sr.? O (o) 

and 

d£?(0) 9S« „(0) 

-£n = =v*'*-« + -^->- < 4 - 28) 

By applying the argument in §4.2 to the above equations we obtain the following behavior of these. 

dk x 1 — a 2 

dS^(O) oi{—vu + avu) 

* ~ ky 9 kx . (4.30) 
dk x 1 — or 
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Then the backflow defined in §2 is derived as follows. 

^ ~ \ (4.31) 

zy, 1 — cr 

Prom these equations it can be seen that the backflow is negative in the case of > and the 
magnitude of this becomes large as the spin fluctuation grows (i.e. a has the larger value). These 
considerations indicate that is smaller than in magnitude and A~ 2 < 47re 2 n/m* follows 
this. 

§5. Results of Numerical Calculation 

In this section we present the results of numerical calculations based on the model and approxi- 
mation presented in the previous section. 

The first Brillouin zone is divided into 128 x 128 mesh and the roughness in the following figures 
of the Fermi surface and the velocity is caused by this mesh structure. The hole doping level is 
taken as 5 = 0.15 in the figure with no indication. 

Firstly (a quarter of) the Fermi surfaces in the first Brillouin zone with the hole doping 5 = 0.15, 
calculationed with SC-SOPT and FLEX approximations are shown in Fig. 3. In this figure hsl 
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Fig. 3. Fermi surfaces with the hole doping 8 = 0.15, calculated by using SC-SOPT and FLEX approximations. 



and hs2 mean the hot spots which are the intersection points between the Fermi surface and the 
magnetic Brillouin zone (the line connecting (jr,0) with (0, ir) in this figure). These are the points 
connected by other points on Fermi surface by the vector Q = (ir, w), and cs means the cold 
spot which is the point (on the Fermi surface) most far from the magnetic Brillouin zone. These 
notations are the same as in some papers on high-T c cupratesS It is also known that the Fermi 
surface is deformed by the correlation effect. From this figure it is noted that one with FLEX is 
deformed much to the magnetic Brillouin zone than the other with SC-SOPT. This is because the 
FLEX approximation includes the effect of the spin fluctuation effect much. 

5. 1 Momentum dependence of the current 

From here the only x-component of the vector is shown in the figure. The y-component at 
(k x ,k y ) is identified with the x-component at (k y ,k x ). j£ in SC-SOPT is shown in Fig. 4(a) with 
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renormalized velocity and bare velocity v^ x . From the argument given in §4.1 the current 
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Fig. 4. The momentum dependence of jjk >fjj and in (a) SC-SOPT and (b) FLEX 



jjjj is not expected to be much different from v*j^ ~ z k v kx anc ^ this ex P ec tation is verified. The 
term 5S^,(0)/3/c x in the renormalized velocity is negligible in the case of weak momentum 
dependence of the irreducible four point vertex. To put it more precisely, replacing by v^, 
like in ref.3, is not verified exactly but only approximately with a condition given above on the 
irreducible four point vertex, in FLEX is shown in Fig.4(b). From this figure it can be seen 
that the current is different so much from the renormalized velocity, unlike in the case of SC-SOPT, 
owing to the strong spin fluctuation. There are two reasons of this difference. One is that is 
reduced by the large negative value of which originates from the large . The other is 
that increases due to the positive value of dT l j e (0)/dk x . The latter point manifests itself in 
v kx > v kx at ^ s ^ - Both of these originate from the large spin fluctuation and in this case it is not 
allowed to approximate j\, x — z k v kx' as in SC-SOPT. The above mentioned difference between 
jjk and also numerically verifies the discussion on the backflow in §4.4. 

5.2 A~ 2 

Based on the above behaviors of jt the doping dependence of A -2 is presented. The result of 
SC-SOPT is shown in Fig.5(a). From this figure it can be seen that A -2 , which is proportional 
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Fig. 5. The integrals over the Fermi surface, J* Fg 2 J* V , | v^X, in (a) SC-SOPT and (b) FLEX. Here X = f kx 
z k v k, and v k x - 



to the case of X = , is approximately described by n/m*, which is proportional to the case 
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of X = ZuVu ■ This is because the momentum dependence of the irreducible four point vertex is 
weak as discussed in §3 and §4.1. It is also noted that the value of the case X = vu decreases as 
the hole doping 5 increases. This is because the case X = vu corresponds to n/m, n means the 
effective carrier and increases owing to the increase of the effective carrier density 1 — 5. (5 = 
means the half- filled case and therefore the effective carrier density is 1.) The decrease of n/m as 
5 increases indicates that n/m* = zn/m (the case of X = ) also decreases because the doping 
dependence of z is rather weaker than the change of n. The above results are considered to explain 
the experimental doping dependence of A -2 in the overdoped region because the perturbation 
scheme is valid due to the weakness of the spin fluctuation in this region. The calculations by 
FLEX is shown in Fig. 5(b). Unlike the case of SC-SOPT, it can be seen that the value and the 
doping dependence of the case X = jt is very different from those of the case X = z^v^. 
This is because the spin fluctuation highly renormalizes the current. The doping dependence of 
the case X = zuvu is rather weak compared with the case of SC-SOPT owing to decreasing of 
the renormalization factor in the low doping. On the other hand the value of the case X = jt 
decreases as the 5 decreases because the spin fluctuation grows in the low doping region. This 
behavior of A -2 (the case of X = jjL ) is considered to explain the experimental doping dependence 
in the optimal and the under doped regions. 

While the physical meaning and a justification of the usage of the perturbation theory on the 
overdoped region and the spin fluctuation theory on the optimal and the underdoped regions is 
discussed in §7, we discuss the following two points. One of these is that by discussing the above 
results conversely, it can be said that the peak at the slightly overdoped region not at the optimal 
doping as experimentaly observed suggest that the spin flucutation begins to grow in this region, 
and therefore this is reflected in other quantities like the one-particle spectrum.0) The other is 
about the renormalization factor z. The calculated z by SC-SOPT and FLEX are not smoothly 
connected. This suggests that the higher order terms in the perturbation expansion are needed to 
calculate the doping level itself of the peak in A -2 as discussed in §7, while the explanation on the 
behaviors of A~ 2 on both sides of the peak is not modified. 

5.3 On the slope of A~ 2 at absolute zero 

In high-T c cuprates it is known that at the low temperature A -2 decreases linearly with T because 
of the line nodes in the superconducting gap. If we assume the isotropic case, the value of A -2 at 
T = is given by, 



m* \ 2 

which comes from the fact that A -2 is linear in j*. On the other hand from eq.(4.22) it can be seen 
that the coefficient of T in A~ 2 at the low temperature is given by 

«3 ! + T > ( 5 - 2 ) 



dT m* V 2 
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which comes from the fact that dA~ 2 /dT is square in j* in this case. Then it is expected that if 
the decreasing of A~ 2 is attributed to small 1 + F 2 /2 as 5 decreases, then the rate of decreasing of 
dA~ 2 /dT is expected to be more rapid than A~ 2 owing to (1 + F 2 /2) 2 . This expectation is denied 
by the experimental results (e.g. refs.25 and 26). Here we clarify that this failure in the explanation 
for the experiments is not caused by the failure of the explanation based on the Fermi liquid theory 
but owning to the use of the isotropic model which is unrealistic as understood from the previous 
sections. 

The slope dA _2 /dT at T = comes from the derivative of Yosida function by the temperature 
dY(k;T)/dT. From the fact that the excitations of the quasiparticles at the low temperature are 
mainly produced at the line nodes (which are equal to the cold spots in our system), dY(k;T)/dT 
has a large value at these points. Therefore it is allowed to consider the function jj^ 2 only at the 
cold spots in dA~ 2 /dT|T=o- On the other hand for the value of A~ 2 at T = the whole value of 
v kx^kx on Fermi surface should be considered in practice. 

With the above consideration in mind we present the momentum dependence of jj^ at hole 
doping 5 = 0.10 and 5 = 0.20 in Fig.6(a). From this figure it can be seen that the doping 
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Fig. 6. (a) The momentum dependence of jj^ at 8 = 0.10 and S — 0.20 in FLEX, (b) The doping dependence of 
at the characteristic points hsl, cs and hs2 on Fermi surface in FLEX. 



dependence of the current around k = hs2 is large while those around k = hsl,cs is weak. 
The doping dependences at these points are shown in Fig.6(b). From these figures it can be seen 
that the decrease of A -2 at T = presented in fig.5(b) is caused by the decrease of around 
the points hs2 while jj^ around the points hsl and cs doesn't contribute so much on the doping 
dependence of this quantity. On the other hand the small doping dependence of around cs is 
reflected in the small doping dependence of dA~ 2 /dT. These behaviors of jj^ , with the anisotropy 
included, are considered to explain the experimental results of A~ 2 and dA~ 2 /dT consistently. 

§6. Relation between Magnetic Field Penetration Depth and Superconducting Fluc- 
tuation 

The high-T c cuprates is known to show the pseudogap phenomena in the underdoped region. 
Some authors are tried to explain this phenomenon by relating this with large value of the magnetic 
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field penetration depth. If The point of their argument is that by writing the phase only model the 
phase stiffness is proportional to the inverse of squared magnetic field penetration depth and then 
the large value of A at absolute zero in the underdoped region means that the phase fluctuation is 
dominant in this region. In this section we investigate whether this assertion is correct. 
The GL model is written as 



H = / dr 



^o|^| 2 + Co 2 |(-iV + ^A)^| 2 ) + ^| 



(6.1) 



Here a := p*(0) (p*(0) is the density of states at Fermi level enhanced by 1/z), eo := (T — T c )/T c , 
£o := ( ; ^r~) ' ^0 is the flux quantum and b : = (nT°)' 2 by ^he microscopic calculation with 
an usual model where the attractive force is operated between renormalized quasiparticles. By 
using if; = | -0 1 e 1 ^* this model is reduced to a phase only model, 

okol 2 1 . / . ^o_» i 



X-\T)=( 2 ^\^A. (6.3) 



H = I dr 

Here we put 

\- 2 (T) = C- 

From this expression with the above parameter inserted it is seen that the A -2 is proportional to 
z. This is in contrast to A~ 2 oc z(l + Ff/2) at absolute zero discussed in the previous sections. On 
the other hand, from eqs.(2.2) and (2.5) the magnetic field penetration depth near the transition 
temperature is given by 

1 f dS k 

tftv(T ~ T c ) J FS 27rV(fc)| 

eo term is derived from 1 — Y(k;T) and we neglected the second term in the right hand side of 

eq.(2.5) due to the smallness of 1— Y(k\ T) near T c . From this equation A~ 2 (T ~ T c ) is proportional 

to the renormalization factor z and this is consistent with the eq.(6.3) not with the expression of 

A~ 2 at T = 0. This suggests that Ff term which is one of the particle-hole four point vertices is not 

effective near T c where the particle-particle four point vertex is remarkable, although this notion is 

much simplified and should be elaborate for the exact description. 

The above consideration indicates that the the smallness of the coefficient in the phase only 

model cannot be identified with the smallness of A -2 at the absolute zero even if the temperature 

dependence eo is excluded. This is originated from the fact that there is a non trivial temperature 

dependence in the vertex correction in the superconducting state. Therefore the small A -2 at 

absolute zero doesn't necessarily indicate that the thermal fluctuation around T c is large because 

the renormalization factor in the phase model is not z(l + Ff/2) but z. Then the the small 

A~ 2 cannot be directly related with large thermal fluctuation in the underdoped region although 

it is considered that the thermal fluctuation is large as experiments indicate, but the origin of 

the large fluctuation is composed of the factors like the strong coupling effect and the quasi two 

dimensionality]!' 
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§7. Summary and Discussion 

In this paper we derived the general expression for the effect of the Umklapp scattering on the 
current and discussed the magnetic field penetration depth, particularly of high-T c cuprates, both 
analytically and numerically. 

In §2 a formula for the current carried by quasiparticles which explicitly expresses the effect of the 
Umklapp scattering is derived. This is basically derived from the usual expression for the current 
written by the renormalized velocity with the backflow. However to see the effect of the Umklapp 
scattering the above two quantities are not basic but this effect is written by the integrals over the 
whole energy scale not only over the Fermi surfaces. This is the basic difference between the current 
in the hydrodynamic region and in the collisionless region, and indicates that the backflow term 
is not necessarily the basic quantity in the lattice system. The advantages of our expression for 
the Umklapp term over the notion of the backflow are the followings. We can show that the upper 
limit for the current is the bare velocity and doesn't exceed this value and that the correlation 
effect is unified in the Umklapp term while by using the backflow notion the renormalized velocity 
also includes the correlation effect. 

In §4 the Umklapp term, mainly of the high-T c cuprates, is analytically discussed. In the highly 
overdoped region of the high-T c cuprates the momentum dependence of the irreducible four point 
vertex is weak though it exists, and by neglecting this small term the Umklapp term is written 
only by the self-energy term. Therefore the current is expressed only by the product of the renor- 
malization factor and the bare velocity and then the A -2 is proportional to n/m*. By this fact the 
reason why a misunderstanding of A~ 2 oc n/m* like ref.3 has not been confronted with a difficulty 
so far is understood. In the optimal and underdoped regions it is known that there exists the 
strong spin fluctuation, and the temperature-linear resistivity and the temperature dependent Hall 
coefficient can be derived by using this feature. Therefore it is reasonable to start with the spin 
fluctuation model to discuss the magnetic field penetration depth too. It is found that the strong 
momentum dependence of the irreducible four point vertex induces the large vertex correction in 
addition to the self-energy term and this correction makes the Umklapp term dependent on the 
points on the Fermi surface. Therefore at some points the current is not written by the product 
of the renormalization factor and the bare velocity unlike the case of the highly overdoped region, 
but reflects the strong Umklapp scattering at these points. These considerations indicate that the 
strong antiferromagnetic spin fluctuation can be identified with the strong Umklapp scattering as 
the effect on the current. The other characteristic results are that the backflow always takes an 
opposite sign to the bare velocity and the partial derivative of the self-energy by the momentum 
takes the same sign as the bare velocity with non-negligible magnitude. The latter point indicates 
that the renormalized velocity is not so simple and supports the above assertion that the dividing 
the current to the renormalized velocity and the backflow is not useful. 
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The main results obtained by numerical calculations in §5 are that the doping dependence of A 
in experiments can be explained by using the perturbation scheme in the overdoped region and 
the spin fluctuation theory in the optimal and underdoped regions, and that the relation between 
A" 2 and dA~ 2 /dT at absolute zero indicated by experiments, which was failed to explain by the 
isotropic model,© is explained by using the anisotropic model based on the spin fluctuation theory. 
As regards with the latter point the way to recover the consistency between these two quantities 
based on the isotropic model is indicated in ref.25 by using the angle resolved photo-emission 
spectroscopy (ARPES) experiments. It is summarized as that ARPES experiments indicate that 
the form of the gap function deforms from the d x i _ y 2 -pairing with under doping and this suggests 
that the magnitude of the gap around the cold spot decreases and then the excitations around 
this point increases with underdoping. Although this behavior may be possible in the underdoped 
region, the calculated gap with the FLEX approximation suggests that the deformation from the 
d x 2_ y 2 -pairing is not so drastic as this experiments suggests. While one of the reasons for this is 
the inaccuracy of the approximation, the other is that the ARPES experiments at the cold spot 
where the dispersion is sharp, cannot be accurate to discuss the form of the gap function around 
the node because of the limited resolution in ARPES. In spite of these facts the anisotropy of the 
system is not obviously negligible and it is needed for considering the qualitative estimation for 
dA- 2 /dT. 

The explanation on the doping dependence of A -2 needs some notions. It is desirable not to 
rely on the specific approximations, but this cannot be done. Then it is needed to consider what 
is physically the most appropriate method. Because we stand upon the Fermi liquid theory, the 
perturbation scheme should be taken at first. In high-T c cuprates the on-site Coulomb interaction 
in the Hubbard model is large, and then it may be considered that all higher order terms is effective 
in the expansion. The reason why the some low order terms is sufficient in the overdoped region and 
the spin fluctuation theory is appropriate (only approximately) in the optimal and the underdoped 
region, is explained as follows. In the study on the single impurity system, it is shown that in 
the perturbation expansion by the on-site Coulomb interaction U, the coefficient of U n (n is the 
order of the expansion) becomes smaller as n becomes large.© This indicate that it is practically 
sufficient to take only the low order terms in U, even if U is somewhat large, in contrast to the 
mean field treatment. The notable point is that there is no momentum dependence in the single 
impurity system. This indicate that it is reasonable to take the low order perturbation scheme as 
an approximation in the overdoped region where the momentum dependence of the vertex is rather 
weak owing to its filling. With underdoping the momentum dependence of the vertex becomes 
strong because the system approaches to the antiferromagnetic phase. Therefore in this case taking 
only the low order terms is not appropriate owing to the fact that the cancellation between higher 
order terms is not effective, but the specific mode is considered to be effective up to higher order 



21 



terms and this leads to the antiferromagnetic spin fluctuation. This is the reason why we take the 
spin fluctuation theory as an approximation in the optimal and underdoped regions. It is also the 
reason why FLEX approximation is inappropriate in the overdoped region because other higher 
order terms which should cancel the RPA type terms in the case of the weak momentum dependence 
are not included. The above qualitative discussions should be verified by explicitly calculating the 
higher order terms with various doping levels. This is an important feature problem. 

The other notions with the FLEX approximation is that this approximation is not good to 
describe the Mott transition® The main fault of this approximation is lack of the Hubbard 
peaks in the density of states© However it is the momentum dependence at the low energy that 
is important for our explanation on the current and the magnetic field penetration depth, and 
therefore main results in this paper is not changed qualitatively by the existence of the high-energy 
structure like the Hubbard peak while the improvement of the density of states would slightly affect 
the quantitative results. 
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Appendix: Expressions at Finite Temperature 

At finite temperature the corresponding equation of eq.(2.20) is (w^(e) := — 1 J| — v k^> 

dGUe') 



1 




'k' 

de' 



i v k' ~ v k) 



+ 2i jk' 



de' 



+ 3, Vde.e' 



de' J l-0£j*(e / )/0e / k,k 



1 ; T^l(e,e^Gi^ ^ 



k y 'l-d^{e')/de r 



Ik' 2T coshV/2T k,k' 
Here the meanings of (11), etc. of 9 are same as in Eliashberg's paper© and these are 



l-d^,(e')/de'_ 
(A-l) 



(e, e') = tanh^[r^,(e, e>) + e')\ + 2fco S h^A«,(e ) e'), (A-2) 



2T k,k 



k,k 



2T k,k 



and 



duo 



\w=0 



,(12) 



2T CO sh 2 e'/2T k,k 
2T coshV/2T [F k,k'^ €>) + F k,k' (e ' e ' )] 



9 (13) (e e' 
^k,k' [ ' 



k,k 



2T k,k 



(A-3) 
(A-4) 



tanh^[r^,(6, e>) + T^\,(e, e>)\ + 2icosh^ A^^e, e'). (A-5) 
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Then 



e e \(°) *v -^ fc fc v ' e j~ e A ( ) ,v ,/\ fc v y fc v 7 



ro th- -A w (V f 'l fc fc roth— — A 1 > (f f')- 



2T fc,fc' v ' ; de' (e')/de' 2T k,k' y ' ; <9e' 1 - <9£« (e')/de'' 

(A-6) 

Here we put r^'?\(e,e') = Rer^'?], (e, e') + iA^'/^e, e') as the four-point vertex which has a dis- 

AC, AC AC, AC AC, AC 

continuity along the cut Im(e — e') = or Im(e + e') = for the case of (a) or (b) in superscript, 
respectively. u^(e) has the following form and it is seen that this quantity vanishes in the absence 
of Umklapp process. 

r f-f' /a OGfAe') f + f > dGfAe') 

- tanh 2T^ [r L a !fc' (e ' e,) +r fc!fc' (e ' e ' )]ImG fc' (e/)}( ^ - ( a - 7 ) 

r r r dx cosh ^ff^ , 2 

~~ Jk'JqJ ~ cosh^cosh^cosh^ ' A;,fc',q e ' W ' x ) I 

ImGf_ a (e-^)ImGf,(e / )ImGf, + „(e , + x) / , A N 

x — ^ — {x -u-i5r {Vk « + v «+« ~ v « ~ Vk) - (A ' 8) 

Prom eq.(A.7) to eq.(A.8) we used the spectral representation 

f dx Imgg( 



G k 



rd_x^ x) 

W J vr x - e - id V ; 



(here (5 is positive infinitesimal quantity) and the trivial relation between hyperbolic functions 

x ,w-eW f' 1 x + e'\ cosh^±^ /k 

cosh— - tanh — — tanh— - tanh — — — = ^ 7— (A- 10) 

2T 2T J \ 2T 2T J cosh^ cosh ^ cosh ^ v > 

Prom eq.(A.6) and eq.(A.8) it is obvious that in the continuum wjj,(e) = also holds at finite 
temperature as derived at absolute zero in §2. By taking the SC-SOPT and FLEX approximations 
as example the validity of the above equation can be confirmed. In the latter approximation it is 
needed to use the exact relation lmVq(uj) = Wq{uj)lm.Xq{io) ■ 
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